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(N ; ABSTRACT 

O 

CD 

Q ■ High resolution spectroscopy of AE Aquarii with the XMM-Newton RGS has 

enabled us to measure the electron number density of the X-ray-emitting hot 
^ ■ plasma to be ~ 10^^ cm~^ by means of intensity ratios of the He-like triplet of 

Nitrogen and Oxygen. Incorporating with the emission measure evaluated by the 
EPIC cameras, we have also found a linear scale of the plasma to be ^ 5 x 10^'' cm. 



> 
0^ 



•/^ ■ Both these values, obtained model-independently, are incompatible with a stan- 

^ ■ dard post-shock accretion column of a magnetized white dwarf, but are naturally 

^ ■ interpreted as the plasma being formed through interaction between an accre- 

O ■ tion flow and the magnetosphere of the white dwarf. Our results provide another 

^ . piece of evidence of the magnetic propeller effect being at work in AE Aqr. 

O \ Subject headings: binaries: close — novae, cataclysmic variables — stars: indi- 

\ vidual (AE Aquarii) — plasmas — X-rays: stars 



1. Introduction 

AE Aqr is a close binary system composed of a magnetized white dwarf rotating at a 
period of 33.08 s (Patterson 1979) and a Roche-lobe filling K3IV secondary star orbiting at a 
period of 9.88 h (Welsh, Horne, & Comer 1993). The masses of the primary and the secondary 
are reported to be Mi = 0.79 ± O.I6M0 and M2 = 0.50 ± O.IOM0 (Casares, Mouchet, 
Martinez-Pais, & Harlaftis 1996), and the orbital semi-major axis length is a = 1.8 x 10^^ cm. 

AE Aqr has been known as one of the most enigmatic magnetic Cataclysmic Variables 
(mCVs) on various aspects, including large optical flares and flickering (Patterson 1979), 
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large radio flares (Bastian, Beasley, & Bookbinder 1996), TeV 7-ray emissions (Meintjes et 
al. 1994) and so on. Although the orbital period of this long predicts existence of an accretion 
disk, an optical spectrum of AE Aqr does not seem as such, but only shows a broad single- 
peaked Ha emission whose centroid velocity was found to lag behind the white dwarf orbit by 
some 70° (Welsh, Horne, & Gomer 1993). Hard X-ray emissions from magnetic CVs (mCVs) 
originate from the post-shock plasma that takes place close to the white dwarf surface in the 
accretion column. Although the plasma temperature is a few tens of keV in general (Ezuka 
& Ishida 1999; Ishida & Fujimoto 1995), that of AE Aqr is measured to be as low as ~ 3 keV 
(Eracleous 1999; Choi, Dotani, & Agrawal 1999). All these results, together with discovery 
of a steady spin-down of the white dwarf at a rate P — 5.64 x 10~^^ s s~^ (de Jager, Meintjes, 
O'Donoghue, & Robinson 1994), lead Wynn, King, & Horne (1997) to draw a picture that 
the accreting matter from the secondary trapped at the magnetosphere does not accrete onto 
the white dwarf but is expelled from the binary owing to magnetic torque (Welsh, Horne, & 
Gomer 1998; Ikhsanov 2001). This is so-called the magnetic propeller effect. 

In this paper, we present new observational results that strongly indicate the magnetic 
propeller being at work in AE Aqr, by means of a model-independent density diagnostics 
with the He-like triplet of Nitrogen and Oxygen fully resolved by the XMM-Newton RGS for 
the first time. 



2. Observations 

AE Aqr was observed with XMM-Newton (Jansen et al. 2001) on 2001 November 7-8. 
The observation log is given in Table 1. The entire observation consists of two parts. During 
the former ~10 ks, both the EPIC pn (Striidcr et al. 2001) and MOS (Turner et al. 2001) 
had been switched off. The RGS (den Herder et al. 2001), on the other hand, had been 
normally operated throughout the observation. As a result, 27 ks data are available for the 
RGS and 17 ks for the EPIC pn/MOS. In Fig. 1 shown are the light curves of pn, MOSl, 
and RGSl. 

For the data reduction, we use the SAS version 6.0.0. We adopt circular apertures of 
54'.'7 and 50'.'6 in radius as the source photon integration region for pn and MOS, respectively. 
The background photons are accumulated from concentric annuli, being set out of the source 
regions, with three times larger outer radii. Source photons are somewhat piled up in the 
pn data especially during the flare. We thus have excluded the central circular region with 
a diameter of 16" from the analysis of pn data. 
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3. Data Analysis 

3.1. The EPIC Spectra 

The time-averaged energy spectra extracted separately from the MOS and pn data 
are shown in Fig. 2. A number of Ka emission hnes of H-hke and/or He-hke Ne, Mg, 
Si, S, and Fe can be recognized. Coexistence of these hnes indicate that the X-ray emission 
originates from a multi-temperature optically thin thermal plasma (Choi, Dotani, & Agrawal 
1999; Eracleous 1999). Wc thus have adopted a multi-temperature VMEKAL model (Mewe, 
Kaastra, & Liedahl 1995) to fit these spectra. We have added a new temperature component 
one by one until the addition of a new component does not improve the fit significantly on 
the basis of F-test. To this end, we have arrived at a four temperature vmekal model 
with common elemental abundances undergoing photoelectric absorption with a common 
hydrogen column density. The best-fit parameters are hsted in Table 2. Emission measures 
of the four continuum components with the temperatures of 0.14, 0.59, 1.4, and 4.6 keV are 
1.3, 3.6, 2.7, and 5.3 x lO^^cm"^, and 1.29 x lO^^cm"^ in total, for the assumed distance of 
100 pc (Welsh, Horne, & Oke 1993). The highest temperature kT — 4.6 keV is considerably 
lower than any other mCVs. Also, the fiuorescent neutral iron Ka emission line at 6.40 keV 
ubiquitous among mCVs (Ezuka & Ishida 1999) is absent. The upper limit of its equivalent 
width is obtained to be < 88 eV. Thanks to a large effective area of XMM-Newton, the 
abundances of the elements from N to Ni are obtained. They are generally sub-solar (Anders 
& Grevesse 1989), except for N which is more than three times the solar value. 



3.2. Line Intensities from RGS Spectra 

In Fig. 3 we have shown the RGS spectra. The Ka emission lines from Nitrogen through 
Silicon in the H-like and/or He-like ionization states can easily be recognized. Figure 4 (a)- 
(c) are blow-up of energy bands of the He-like triplets of N, O, and Ne. In the next column 
shown are the spectra predicted by the best- fit four temperature vmekal model (table 2) , 
being convolved by the RGSl 1st order response function. Note that the vmekal model 
represents emission spectrum from a thermal plasma in the low density coronal limit. One 
can easily see, for N and O, that the intensity of the intercombination increases from the 
coronal limit by consumption of the forbidden. This behavior is interpreted as a high density 
effect; if the electron density exceeds a certain critical value inherent in each element, one 
of the two electrons excited to the upper level of the forbidden line ('^Si) is further pumped 
by another impact of a free electron up to the higher level ^P2, i, and is then relaxed by 
radiating the intercombination line. The relative intensities of the intercombination and 
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forbidden lines can therefore be utilized as a density diagnostics (Gabriel & Jordan 1969; 
Pradhan & ShuU 1981). 

In order to evaluate the electron number density Ue, we begin with evaluating intensities 
of the He-like triplets. For this, we utihze the four temperature vmekal model that provides 
the best-fit to the EPIC pn and MOS spectra. We have fixed the hydrogen column density 
and the four temperature at the values in table 2 and the relative normalizations of the four 
continuum components. The abundances are also fixed at the best-fit values, except for an 
element to be used as a density diagnostic, for which the abundance is set equal to zero, and 
instead, four Gaussians are added, representing the He-like triplet and the Lja line. The 
line intensities thus obtained for N, O, and Ne are summarized in table 3, and the best-fit 
results are displayed in Fig. 4 (a)-(c) as the histograms. The ionization temperature kTi 
calculated from the intensity ratio between the Lya and r are obtained to be 0.16, 0.30, and 
0.34 keV for N, O, and Ne, respectively. 

Given the line intensities of the triplets, we have carried out density diagnostics by 
means of the intensity ratio f /{r + i). In Fig. 4 (g)-(i) shown are theoretical curves of the 
ratio f /{r + i) versus the electron density, which is calculated by means of the plasma code 
SPEX (Kaastra, Mewe, & Nieuwenhuijzen 1996) at the ionization temperature kT^ of each 
element. In each panel we have also drawn a range of the intensity ratio allowed from the 
data and the resultant density range as a box. The electron densities are obtained to be 
1.4 X IQio - 1.3 X 10"cm-2, 4.0 x 10^° - 6.8 x lO^^cm-^ and < 9.3 x 1012^^-2 for N, O, 
and Ne, respectively. As a crude approximation, the electron number density of the plasma 
is fie ~ lO^'^cm^. 

4. Discussion and Conclusion 

It is of great importance to note that the resultant density ~ lO^^cm^ is smaller 
by several orders of magnitude than the conventional estimate in the post-shock accretion 
column of mCVs: — lO^^cm"^ (Frank, King, & Raine 2002). Moreover, combined with 
the emission measure EM = 1.3 x lO^^cm"^ obtained with the EPIC pn/MOS (§ 3.1), the 
hnear scale of the plasma is evaluated to be ip — {EM/nlY^^ ~ 5 x 10-'^'^ cm, which is much 
larger than a radius of a typical white dwarf and rather close to the orbital scale. We thus 
conclude that the optically thin hot plasma in AE Aqr does not accrete onto the white dwarf 
but rather spreads over the orbit of the binary. Note that this conclusion is derived simply 
on the basis of widely approved atomic physics, and hence, is free from any specific model. 

It is known, however, that photo-excitation due to UV radiation can also pump ^Si 
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electron up to ^P2,i, thereby affecting tfie density diagnostics. Altfiougfi tfiere is no evidence 
of accretion disk in AE Aqr (Welsh, Horne, & Gomer 1993) and the secondary star is a late 
type K3 subgiant, the white dwarf can be a source of the photo-excitation, the emission 
spectrum of which is reported to be a blackbody with a temperature of Tj. = 26, 000 K 
(Welsh 1999). The photo-excitation rate of such a white dwarf from the initial state i = ^Si 
to the final state j = ^P2, i is 

r«(r,) = ^/...i^.^M (1) 

(Porquet et al. 2001; Mauche 2002) where is the effective oscillator strength of photo- 
excitation which is 0.03574 (Nahar & Pradhan 1999), and is a geometrical dilution factor 
to be taken here as {Rwd/^pY where i?^,^ is the radius of the white dwarf. On the other 
hand, the coUisional excitation rate to be compared with P^ is neQij, where, 

%(^e) = ^^^^^e-^^^/^-«7.(Te) [cm^s-] (2) 

is the rate coefficient of the electron impact excitation (Mewe & Schrijver 1978). In this 
equation, uji{= 3) is the statistical weight of the lower level, Eij is the energy difference be- 
tween the two levels, and jij is the Maxwellian-averaged collision strength which is tabulated 
as a function of Tg in Pradhan, Norcross, & Hummer (1981). Adopting Tj listed in table 3 
as Te, we compared eqs. (1) and (2) and obtained P^j ~ 0.07 UeQij for both N and O. Hence, 
the photo-excitation effect can be neglected in the first order approximation. 

We remark that the density obtained in § 3.2 should be treated as an upper limit, in the 
case that the UV radiation is stronger than the estimation above. Even if so, the conclusion 
above need not to be changed because a lower density, and hence, a larger geometrical scale 
are required for the plasma. 

The low density and the large scale of the plasma naturally lead us to invoke the magnetic 
propeller model, in which blobby accreting matter originally following a ballistic trajectory 
from the inner Lagrangian point is gradually penetrated by the magnetic field of the white 
dwarf, and is finally blown out of the binary due to interaction with the rapid rotating 
white dwarf magnetosphere (Wynn, King, & Hornc 1997). The magnetic propeller model is 
advantageous in explaining various characteristics of AE Aqr, such as the spin down of the 
white dwarf, the velocity modulation of Ha line. The low plasma temperature compared 
with other mCVs can be attributed to a halfway release of the gravitational potential energy 
of accreting matter down to the magnetosphere. We befieve that the density diagnostics 
presented here add another piece of evidence to support the magnetic propeller effect being 
at work in AE Aqr. 
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Welsh, Horne, & Gomer (1998) estimated that the velocity of the out-flowing plasma 
can be as high as several hundred km at most. If so, a Doppler shift of the N and O lines 
may be detected through phase-resolved spectral analysis. We thus have made RGS spectra 
of the first and second flares separately, during the time interval of 3,000-8,500 s and 19,000- 
26,000 s in Fig. 1, respectively, and evaluated the central energy of Oxygen Lja. The result 
is negative with the line central energies of 653.4 ± 0.5 eV and 653.4 ± 0.4 eV, respectively, 
which are fully consistent with the laboratory value. We have further spht the data of 
the second flare into evenly segregated three pieces, but the central energy distributes in the 
range 653.2-653.4 eV. Note, however, that we have obtained a finite value for a 1-a line width 
of 2.4 ± 0.6 eV and 1.7 ± 0.6 eV for the first and second flares, respectively, which are larger 
than the natiual width (~ 0.001 cV) or the thermal velocity of Oxygen ion (~ 70 km s~^ 
or ~ 0.2 cV estimated from fcTj = 0.3 kcV). Since the ionization temperature of Oxygen is 
lower than the highest temperature of the plasma (table 2) by an order of magnitude, and 
also from the fact ip 5 x 10^° cm, the Oxygen Lya line emanates from a region far out of 
the magnctosphere where the velocity collimation is already dissolved. 

In order to finally confirm the magnetic propeller effect being at work in AE Aqr, it is 
important to detect a bulk velocity shift of an emission line expected from the plasma flow 
in the vicinity of the magnctosphere. Since the highest temperature of ~ 5 keV, we expect 
this can be done by the iron Ka hne with the Astro-E2 XRS. 
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Table 1. Observation Log 



Obs. ID 


Instrument 


Data Mode 


Obs. Start (UT) 


Obs. End (UT) 


Exp.(s) 


0111180601 


RGSl 


Spectroscopy 


2001-11-07 20:06:35 


2001-11-07 22:53:33 


10018 




RGS2 


Spectroscopy 


2001-11-07 20:06:35 


2001-11-07 22:53:28 


10013 


0111180201 


MOSl 


Large Window 


2001-11-07 23:06:36 


2001-11-08 03:42:09 


16533 




M0S2 


Large Window 


2001-11-07 23:06:35 


2001-11-08 03:42:08 


16533 




pn 


Full Frame 


2001-11-07 23:45:53 


2001-11-08 03:38:17 


13910 




RGSl 


Spectroscopy 


2001-11-07 23:00:19 


2001-11-08 03:45:39 


17120 




RGS2 


Spectroscopy 


2001-11-07 23:00:19 


2001-11-08 03:45:36 


17117 
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Table 2. The best-fit parameters of the 4 temperature vmekal model being fit to EPIC pn 

and MOS spectra 



Parameter 


Value 


Element 


Abundance^ 


iVn (1020cm-2) 


0.0»_^ 20 


N 


q cri +0.92 
'-'•'-'-'--0.81 


kTi (keV) 


4.60t°i° 





74+0-i'^ 

'"±-0.23 


m (keV) 


-1 9-1+0.13 

-"-■"^-■--o.os 


Ne 


0.43l°:i 


kTs (keV) 


0.59lH^ 


Mg 


0.701°:}^ 


kT^ (keV) 


0.141^^ 


Si 


o.8ii°:i^ 


AT.b (10-3) 


4 41^+0-41 

^•^^-0.44 


S 


A 70+0.20 

'-'■'"J-0.18 


iV2'^(10-3) 


9 9C-+0.52 
^•^'-'-0.51 


Ar 


0.21 (< 0.89) 


iV3'^(10-3) 


3.04l°ll 


Ca 


0.19 (< 1.11) 


iV4'^ (10-3) 




Fe 


-0.06 


iVpn/iVMOS^ 


1.081°:°! 


Ni 


-r 97+0.57 
-■-•^ ' -0.50 


(d.o.f.) 


1.22 (992) 







'^Solar Abundances (Anders & Grevesse 1989) 

'^Normalization of the VMEKAL component obtained 
with pn camera in a unit of 10^^'^/47rD^ J UeUH dV, 
where D [cm] is the distance to the target star. 

■^Ratio of continuum normalizations. 

Note. — All the errors are at the 90 % confidence level. 
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Table 3. Intensities of Lya and He-like triplet of Nitrogen, Oxygen and Neon with RGS 



Nitrogen Oxygen Neon 





Energy^ 


Norm'^ 


Energy*^ 


Norm'^ 


Energy^ 


Norm'' 


Lya 
r 

i 
f 


0.50032 
0.43065 
0.42621 

0.41986 


1.85;°il 

0.97l°:i? 

"ST+o-i^ 

'-'■<3'_0.17 


0.65348 
0.57395 
0.56874 
0.56101 


o co+0.34 
O.UO_Q 35 

0.851^:^^ 
0.37l^:i^ 


1.0215 
0.92195 
0.91481 
0.90499 


o.9oi°:i 

0.65l°i^ 

'-'•'-'^-0.18 


xl (d.o.f.) 


0.16 ±0.02 
1.36 (115) 


1.14 (69) 


U-<J^-0.07 

1.21 (63) 



'^ine central energy in keV. 

''ine normalization in a unit of 10~^ photons cm~^ s~^. 

'^Ionization temperature in a unit of keV, evaluated by the intensity ratio 
between Lya and r. 

Note. — All the errors are at the 90 % confidence level. 
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AE Aquarii Light Curves (Background not subtracted) 




Time [sec] 



Fig. 1. — The light curves of EPIC pn and MOSl in the band 0.2-15 keV, and RGSl in the 
band 5-35 Awith a bin size of 256 sec. The source integration region is a circle with a radius 
of 57'.'4 and 50'.'6 for pn and MOSl, respectively. The RGSl light curve is created with the 
1st order photons. 




Fig. 2. — The averaged spectra of EPIC pn (green), MOSl (black) and M0S2 (red) overlaid 
with the best-fit four temperature vmekal model as the histograms. The model components 
are drawn with the dotted lines for MOSl and M0S2. See table 2 for the best-fit parameters. 
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RGSl & 2 Time— averaged spectrum 




10 30 30 

Wavelength [angstrom] 



Fig. 3. — The RGS spectra averaged over the entire 27 ks exposure. The black and red points 
represent the data from RGSl and RGS2, respectively. Identifications of the H-line and He- 
like Ka lines of N through Si are shown with dotted-dashed line. The other unidentified 
lines are mainly those associated with Fe-L transitions. 
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Observed Coronal model Denstiy Diagnostics 




0.9 0.92 0.94 0.9 0.92 0.94 10 15 

Energy [keV] Energy [keV] log [cm~^] 



Fig. 4. — Density diagnostics by means of the He-like triplet. Panels of observed and coronal 
limit spectra around the He-like triplet, and the density diagnostics are arranged from left to 
right columns. The cases of Nitrogen, Oxygen, and Neon are arranged from upper to lower 
raws. Black and red data points and histograms are from RGSl and RGS2. The relative 
intensities of the intercombination (marked with 'i') and forbidden (marked with 'f ') lines are 
undoubtedly inverted for N and O. In the panels (g)-(i), comparisons of theoretical curves of 
the intensity ratio f/{r + i) versus the electron density Ue with the observed intensity ratio. 
In drawing the theoretical curves we have used the plasma code SPEC (Kaastra, Mewe, & 
Nieuwenhuijzen 1996) at the ionization temperatures kTi of 0.15, 0.30, and 0.34 keV for N, 
O, and Ne obtained from the intensity ratio between Lya and r. 



